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THE LINEAR DEPENDENCE OF FUNCTIONS OF SEVERAL VARIABLES, 

AND COMPLETELY INTEGRABLE SYSTEMS OF HOMOGENEOUS 

LINEAR PARTIAL DIFFERENTIAL EQUATIONS* 

BT 

GABRIEL M. GREEN 

Inteodtjction 

In the case of n functions of a single variable, the vanishing of the wronskian 
is the most familiar criterion for their linear dependence. The wronskian 
also plays an important part in connection with the theory of a single ordinary 
homogeneous linear differential equation of the nth. order, even in questions 
which do not concern the linear dependence of solutions of the equation. 
It is the purpose of the present paper to generalize the fundamental facts 
connected with these topics, to the case of functions of several variables. 

The characteristic property of an ordinary homogeneous linear differential 
equation is that any solution is expressible linearly, with constant coef- 
ficients, in terms of a fundamental set of solutions. The natural generalization 
to several independent variables is afforded by the completely integrable 
system of homogeneous linear partial differential equations, in one dependent 
variable, any solution of which is likewise linearly dependent on a finite number 
of solutions of the system. It is such systems which are discussed in the 
latter part of this paper; the subject of the linear dependence of functions of 
several variables is developed in the first part, chiefly with a view to its sub- 
sequent application to completely integrable systems. 

In our discussion of linear dependence, we consider throughout n functions 
2/1 J 2/2 . • • • > Vn oi p independent variables ui, U2, • • • , Up. Functions and 
variables may be either real or complex; if any of the variables be complex, 
we suppose the functions to be analytic in those variables. All of the main 
theorems, however, are stated and proved under the supposition that the 
independent variables are real; the modifications — usually simplifications and 
omissions — are easily supplied for cases in which one or more of the inde- 
pendent variables is complex. The n functions, then, of the real variables 

* Presented to the Society in two communications under different titles, October 31, 1914, 
and October 30, 1915. The results concerning the linear dependence of functions of-several 
variables were read on the latter date. 
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Ui, Ui, • • • , Up, exist, together with such of their derivatives as come into 
the discussion, in a jJ-dimensional region, open or closed, which for simplicity 
may be taken to be a p-dimensional parallelopiped. This region we denote 
throughout by A . We shall suppose that all partial differentiations taken 
with respect to more than one of the independent variables are independent 
of the order in which the differentiations are effected. Thus, by the existence 
of d^ y/dudv we imply the existence of d^ y/dvdu and the identity of. these 
two derivatives. 

In the case of functions of several variables, there is no one determinant 
which may properly be taken as the generalization of the wronskian. In 
fact, since there is no relation between the number of functions and the number 
of independent variables, no definite form can be prescribed for such a deter- 
minant; and moreover, the difficulty in this direction is increased by the 
fact that we do not assume our functions to be analytic. However, in the 
general criteria which we obtain concerning the linear dependence of n func- 
tions, a matrix takes the place of the wronskian. This matrix contains n 
columns, and each row is formed of the same derivative of the n functions. 
Evidently, the identical vanishing of all the n-rowed determinants of such a 
matrix is a necessary condition for the linear dependence of the n functions. 
The sufficient conditions at which we arrive are expressed in terms of the 
vanishing of the n-rowed determinants of a matrix of the kind just described, 
and are therefore also necessary — a point which we shall not mention in the 
body of the paper. The case of n functions of a single variable is of course 
included in our general theory; the matrix for such functions may have n 
or more rows, and if it has only n the determinant of the matrix is precisely 
the wronskian of the n functions. 

So far as we are aware, the only completely integrable systems which have 
been studied are of two kinds. Of these, the only non-analytic ones are of 
the first order, and are in fact simply systems of total differential equations.* 
For systems consisting of equations of order higher than the first, an analytic 
character has been presupposed t in establishing the existence of solutions. 
By a mere change of notation, however, we cause such a system to take on 
the form of a system of total differential equations, so that the existence of 
solutions may be proved in the non-analytic case for completely integrable 
systems consisting of equations of unrestricted orders. 

Although for a set of functions of several variables there is in general no 
distinctive determinant, formed from partial derivatives of the functions, 

* Cf., for instance, Mayer, Ueber unbeschrankt integrable Systeme von linearen lolalen Differ- 
entialgleiehungen und die simultane Integration linearer partieller LHfferentialgleichungen, 
Mathematische Annalen, vol. 6 (1872), pp. 448-470. 

t Cf . C. Riquier, Les systemes d'equationa aux derivees partielles, Paris, Gauthier-Villars, 
1910. 
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which might play the role of a wronskian, we have been able to define a wron- 
skian for a fundamental set of solutions of a completely integrable system. 
By means of this wronskian, if we may so term it, we find certain properties 
of homogeneous linear completely integrable systems quite analogous to 
those of an ordinary differential equation, linear and homogeneous, of the 
nth order. An important result in this connection is a proof that a large 
and important class of completely integrable systems of the kind considered 
may be put into a so-called canonical form by a certain transformation of the 
dependent variable. The theorem is analogous to the familiar one in the 
case of an ordinary differential equation 

according to which the equation can be reduced in essentially but one way to 
a unique form for which the coefficient pi is zero, by a transformation of the 
dependent variable alone of the form y = '\{x)y. 

The reduction, just referred to, of a completely integrable system to a 
canonical form, is very useful in certain geometric questions. Completely 
integrable systems have been important in differential geometry, but not 
until recently has systematic use been made of them in connection with that 
branch of the subject for which they are most peculiarly fitted, viz., projective 
differential geometry. We owe to Wilczynski a very general method for 
dealing with questions in that field; in fact, according to this method, the 
projective differential geometry of a configuration is equivalent to the theory 
of the invariants and covariants, under certain continuous groups, of a com- 
pletely integrable system, the configuration being defined by a fundamental 
set of solutions of the system. The continuous groups which give rise to the 
invariants and covariants are generally broken up into two parts: a trans- 
formation of the independent variables, and a transformation of the de- 
pendent variable. The coefficients of the canonical form are invariant under 
this latter transformation, and in terms of these coefficients and their deriva- 
tives may be expressed all quantities having the same invariantive property. 
It is therefore highly desirable to have a method for computing these semin- 
variants, as they are called; we show that this can be done in a purely mechanical 
way for certain completely integrable systems. In the various particular prob- 
lems studied by Wilczynski and others, this method is the one which has 
always been followed. 

As an illustration, we show in § 7 how the foregoing may be applied to the 
theory of curvilinear coordinates in space of n dimensions. 



Trans. Am. Math. Soc. 38 



486 



G. M. GREEN: 



[October 



1. The fundamental theoeem 

Before considering the linear dependence of n functions, it will be con- 
venient to prove the following theorem for two functions: 

Theoeem I. Let yi and y^ be functions of the p independent variables ui, U2, 
• • • , Up for which all partial derivatives of the first order, dyi/duk, dyi/duk, 
{k = 1,2, • • • , p) exist throughout the region A . Suppose, further, that one 
of the functions, say yi, vanishes at no point of A . Then if all the two-rowed 
determinants in the matrix 



yi, 


2/2 


dyi 


dy2 


dUi' 


dui 


dyi 

dU2' 


dy2 

duo 


dyi 


dyi 


dup' 


dUp 



vanish identically in A , yi and y^ are linearly dependent in A , and in fact 



For, if 



we have on dividing by y\ 



whence 



^2 = cyi . 



dyi 
y^du,- 


^2/1 n 


duk\ 


yj 


2/1 


yi = cyi, 



(A; = l,2, 



(fc = l,2, 



p), 



P), 



where c is a constant. 

We might of course have made the statement of the hypothesis somewhat 
weaker, at the sacrifice of smoothness of statement, since all the two-rowed 
determinants of the matrix vanish if all those vanish in which the first row of 
the matrix constitutes a row. 

We note at once that in the case of functions of more than one variable 
there does not exist a single determinant which may properly be called the 
wronskian. We might of course speak of the two-rowed determinants of the 
above matrix as wronskians, but in the case of more than two dependent 
variables, which we shall now consider, new difficulties arise. First of all, 
there is no connection between the number of functions and the number of 
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independent variables. Moreover, the various partial derivatives of the 
functions may exist in a very unsymmetric way, that is, a function may for 
instance be analytic in some of the variables, and may have only a few deriva- 
tives in the others. Consequently, there is not even a set of determinants, 
the structure of which we may describe once for all, which might take the 
place of the single wronskian in the case of functions of one variable. As a 
matter of fact, we shall see presently that very useful sufficient conditions 
for linear dependence may be stated, in which the partial derivatives with 
respect to some of the variables need exist only to the first order, while those 
with respect to other variables may be of higher orders. 

We proceed to the fundamental theorem, of which Theorem I is a special 
case. Let yi, y^, • • • , j/n be functions of the p independent variables u-i,Ui, 
• • • , Up. Let us denote by y^P , y^P , etc., partial derivatives of yi , of any 
order or kind whatever. As a matter of fact, it will not be necessary to 
specify just what derivative of yi is denoted by the symbol y'-p . If each one 
of the set of functions yi {i = 1 , 2 , • • • , n) possesses a certain partial deriva- 
tive, we shall say that the set of functions possesses that derivative. Thus, if 
d'^ yi/dui dui , (i = 1 ,2, • • • , n) exist, the set of functions possesses the 
derivative d^ y/dui du^ . 

We shall frequently have occasion to consider matrices of the form 



2/1. 


Vt, 


••, yr 


y\'\ 


y'P, 


■■. l/p 


yT, 


yT, 


••, y?' 


y\'\ 


iH\ 


••, y'r' 



in which the first row consists of the functions 2/1,2/2, • • • , yr, and the other s 
rows of derivatives of these functions. We shall denote throughout the 
matrix just written by the symbol 

M,{yi, 2/2, •••, yr). 

We may now state the fundamental theorem: 

Theorem II. Let the set of n functions 2/1,2/2, • • • , yn of the p independent 
variables Ui, u^, • • • , Up possess enough partial derivatives, of any orders what- 
ever, to form a matrix 

M = Mn-2(yi, 2/2, ••• , yn), 

of n columns and n — 1 rows, in which at least one of the (n — 1 )-rowed deter- 
minants, say 
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, 2/n-l 


iR\ y'P, ■■■ 


, 2/i'li 


^»-2)^ ^.-2)^ ... 


, if:-? 
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Wr, 



vanishes nowhere in A . Suppose, further, that all of the first derivatives of 
each of the elements of the above matrix M exist, and adjoin to the matrix M such 
of these derivatives as do not already appear in M , to form the new matrix 

M' = M,(yuy2, •••,yn), 

which has n columns and at least n rows, so that q ^ n — 1 . Then if all the 
n-rowed determinants of the matrix M' in which the determinant Wn is a first 
minor vanish identically in A , the functions 2/i .' 2/2 > • • • > 2/n ^J'"^ linearly dependent 
in A , and in fact 

yn = Ciyi+C2y2+ ••• + Cn-l ?/n-l , 

the c's being constants. 

Let us denote by Wi, W^, ■ ■ • , Wn the (n — 1) -rowed determinants of 
the matrix M , in such wise that Wi is the determinant which does not con- 
tain the ith. column of M . If we write yT = yt, we have by the hypothesis 

(1) w,yf + W,^p+---+W^yf^O (j = O.l, ... ,q). 

1 of these equations with respect to Uk , 



Differentiating each of the first n 
we obtain the equations 



(2) 



duk ^' ^ duk y- 



U) 



+ 



^ dUk ^" 






difi' 



+ ^1.-^+^2,^^ 







(i = 0,l, ■■■,n-2). 

But since all of the first derivatives of the elements of the matrix M occur in 
the matrix M' , the last part of each of the relations (2) vanishes because 
of (1), so that 

dWi „, , dW^ ,,, , , dWr. 



duj 



.yip 



dui 



■2/f + 



duk 



£' 



(i = o,l, 



n -2). 



Multiplying the first of these by the first minor of y'l" ( = ?/i ) in Wn , the 
second by the first minor of y^x' in Wn, etc., the last by the first minor of 
2/?""^^ in Wn , and adding the results, we find that 



dWi 

duk 



Wn 



dWn 

dUk 



Wi ^0. 
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But this holds for every k; hence by Theorem I, since W„ 4= anywhere in A , 
we have 

In hke manner, 



Wr^l = — Cn-l W„ , 

the c's being constants. Consequently, from the identity 
we have 

Wn{- CiPi - C2y2 - ••■ - Cn-\ ?/n-l + 2/n ) =0. 

Since Wn is nowhere zero in ^ , it follows that 

2/n = Ci 2/1 + C2 2/2 + • • • + C^l 2/n-l , 

whence our theorem.* 

It should be noted that if, as stated in the hypothesis, all the n-rowed 
determinants of M' of which the determinant Wn is a first minor vanish, then 
all the n-rowed determinants of M' vanish. This follows from the non-vanishing 

oiWnA 

As an example of the application of Theorem II, consider the case of four 
functions of two independent variables, say 



Then if in the matrix 



yi =fiiu,v) 



(i = l,2,3,4). 



2/1. 


2/2, 


2/3, 


2/4 


dyi 
du' 


dy2 
du' 


dyi 
du' 


dy^ 

du 


dyi 

dv ' 


dy2 

dv' 


dya 

dv ' 


dyt 
dv 



at least one of the three-rowed determinants vanishes nowhere in the region 
A, a necessary and sufficient condition for the linear dependence of 2/1. yz, 
yz, 2/4 is the identical vanishing in A of all the four-rowed determinants of 

* This proof is a direct generalization of the corresponding one given by Frobenius for 
functions of a single variable, in the case of the wronskian. Cf. B6cher, Certain cases in 
which the vanishing of the wronskian is a sufficient condition for linear dependence, these Trans- 
actions, vol. 2 (1901), pp. 139-149, § 1. 

t Cf. B6cher, Introduction to Higher Algebra, p. 54, for the general theorem on matrices 
involved here. 
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the matrix formed by adjoining to the above one all of the derivatives of 
the second order of 2/1,2/2,2/3,2/4- If, however, at least one of the three-rowed 
determinants in the matrix 



du' 



2/2, 

dyi 
du' 



2/3, 

3y_3, 
du' 



2/4 

^2/4 
du 



5^2/1 ^^2/2 ^^2/3 3^2/4 



du"" ' du" 



dy? 



du" 



vanishes nowhere in A, then a necessary and sufficient condition for the 
linear dependence of 2/1 , 2/2 , 2/3 , 2/4 is the identical vanishing of all the four- 
rowed determinants of the matrix 



2/1. 



2/4, 



92/1 


dyi 


d'yi 


d'yi 


d'yi 


du' 


dv ' 


du" ' 


dudv' 


aw' 


dyi 


dyi 


d'y. 


d'Vi 


d'y4 


du' 


dv ' 


du" ' 


dudv' 


du' ' 



d'yi 



du^ dv 



du^ dv 



so that the existence of even d^ yi/dv^ is not presupposed. 

In the fundamental theorem concerning the linear dependence of n func- 
tions of a single variable, the identical vanishing of the wronskian is shown 
to be sufficient for linear dependence provided the wronskian of n — 1 of the 
functions is not zero anywhere in the interval of the independent variable 
under consideration. This evidently comes as a special case under Theorem 
II. But Theorem II also affords other sufficient criteria in the case of a 
single independent variable. An example will show the general nature of 
these criteria. Let y , z,whe functions of the real variable x in an interval I . 
Denoting by accents differentiation with respect to x, let us suppose that 
the determinant 

y> '< 
y'", z' 

vanishes nowhere in I . Then a sufficient condition for the linear dependence 
oi y , z, w is the identical vanishing in I of the two determinants 



y, 



y 



y > 

.III 



z 



z, 
z', 

III 



w 
w' 

III 



w 



y> 

,111 



z, 
z'", 
z'\ 



w 
w 



w 

III 



.rv 
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However, in any case such as the present one, a very general theorem con- 
cerning linear dependence, due to D. R. Curtiss, is applicable.* 

Before leaving the fundamental theorem, it may be well to consider the 
case of analytic functions. For n analytic functions of a single variable, the 
identical vanishing of the wronskian is a sufficient condition for linear de- 
pendence, and nothing need be demanded as to the non-vanishing of the 
wronskian of w — 1 of the functions. This follows from the fact that for 
analytic functions linear dependence in any subregion implies linear dependence 
throughout their common region of definition. If in Theorem II the func- 
tions involved are analytic, an obvious weakening of the hypotheses may 
be made; for instance, it is sufficient to know that all the n-rowed determinants 
of the matrix M' which contain the elements of Wn vanish throughout an n-dimen- 
sional sub-region of A which contains a single point at which Wn is different 
from zero. Such would be the most general theorem that may be stated in 
this case; however, by strengthening the hypotheses sufficiently we may 
obtain a theorem in which no demand is made concerning the non-vanishing 
of Wn . A very elegant theorem of this kind is the following, which is due 
to O. D. Kellogg:t 

Let yi, 1/2, • • • , ynbe analytic functions of the p independent variables ui, u^, 

• • • , Up — real or complex — having a common region of definition A , and let 
Mq (yi , y2 , • • • , y„) be the matrix formed from the functions yi> y2, • • • , yn 
and all their partial derivatives up to and including those of order n — 1 . Then 
the functions yi, y2, • • • , yn are linearly dependent in A if and only if all of the 
n-rowed determinants of the matrix Mq vanish throughout A . 

2. Two THEOREMS CONCERNING THE LINEAR DEPENDENCE OF U FUNCTIONS 

In the present section we shall prove two theorems, which are generaliza- 
tions of theorems given by BocherJ for functions of a single variable. They 
will be useful later, in connection with completely integrable systems of 
partial differential equations, just as Bocher's theorems are in the case of 
ordinary linear differential equations, of which completely integrable systems 
are a generalization to several independent variables. 

Theorem III. Let 2/1,2/2, • • • , yn be n functions of the variables tii , M2 , 

• • • , Up. Suppose that the set of functions possesses the n — 1 partial deriva- 
tives y'-P , yf-' , • ■ • , 2/""^' , {i = 1 , 2 , ■ ■ ■ , n) , so that the function 

* Cf. Theorem VI, p. 292, of his article, The vanishing of the wronskian and the problem of 
linear dependence, Mathematische Annalen, vol. 65 (1908), pp. 282-298. 

t Nomograms with points in alignment, Zeitschrift ftir Mathematik und 
P h y s i k , vol. 63 (1914), pp. 160-161. The statement there given is equivalent to ours, 
if the proper interpretation is made of Prof. Kellogg's definition of the rank of a matrix. 

% Cf. Theorems V and VI of the paper already cited, these Transactions, vol. 2 
(1901), pp. 139-149. 
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2/ = Ci 2/1 + C2 2/2 + • • • + c„ 2/„ , 
where the c's are constants, has the same derivatives, y'-^^ , y' 



(2) 



„(n-l) 



2/1 » 


2/2, 


••, 2/n 


2/i", 


2/^', 


••, ^^ 


fr'\ 


2/2"-", 


•••, 2/r^' 



pose also that no function of the form of y , other than zero, vanishes together with 
its derivatives y'-^^ , y^^^ , • • • , y'-"'~^^ at any point of A . Then if the determinant 



W 



vanishes at any point {a-i, a^, • • ■ , ap) of A , the functions 2/1,2/2, • • • , 2/n c-^^ 
linearly dependent. 

For, since W {ai, a2, • • • , ap) = , there exist n constants Ci , C2 , • • • , c„ 
not all zero, such that 

Cil/iHai, ••• ,ap) + c^yfiai, •■• , ap) -\- ••• + Cnyn'{ai, • • ■ , ap) =0 

(i = 0,1, •••,n-l). 

so that the function 2/ = ^i 2/1 + C2 2/2 + • • • + c„ 2/» vanishes together with 
its derivatives 2/^" , 2/'^^ , • • • , 2/'""" at a point of A , and must therefore be 
identically zero. 

Theobem IV. Let yi, y^, • • • , yn be a set of n functions of the p variables 
Ui,U2, • • • , Up which throughout A possesses r partial derivatives (r > n — 1) , 



yf, iV, 



y'V 



(i = 1,2, •••,n). 



M = Mm{yi,yi, ■ 


■■ ,yn) 


(m^n - 2), 


M' = M^iyuyi, • 


■• ,yn) 


(3 > m) 



Suppose that from among these derivatives may be selected at least n — 2, say 
tfP > l/P, • • • , 2/r' (m ^ n — 2) , such that all of their first derivatives occur 
among the set yP , i/P , • • • , yV , and are bounded in A . Let 



and suppose that 



is the matrix formed from M by adjoining to the elements of M those of their 
first derivatives which do not already appear in M . Assume also that if for any v 
columns of M {v = n) all the v-rowed determinants containing elements of the 
first row of M vanish identically, then all the v-rowed determinants formed from 
the same v columns of M' , and containing elements of the first row of M , and of 
V — 2 other rows of M , vanish identically* Suppose, further, that no function, 
other than zero, of the form 

* Some remarks concerning this paxt of the hypothesis will be found after the proof of the 
theorem. 
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2/ = Ci 2/1 + C2 2/2 + • • • + CnPn, 

where the e's are constants, vanishes together with its derivatives y^^^ , y^^^ , • • • , y^"^ 
at any -point of A . Then the identical vanishing of all those n-rowed determinants 
of M' in which appear elements of the first row and n — 2 other rows of M is a 
sufficient condition that the functions 2/1 > 2/2 , • • • , ynbe linearly dependent. 

Suppose first that one of the (n — 1 ) -rowed determinants of the matrix M 
containing elements of the first row of M does not vanish identically. We 
may without loss of generality take for this determinant the one which appears 
in the upper left-hand corner of ilf . As in § 1, we denote this determinant 
by Wn- Then at some point (ai, a2, ■ ■ • , ap) oi A , Wn is different from 
zero. But Wn is continuous in A , because all its partial derivatives of the 
first order exist and are bounded in ^*; consequently, it differs from zero in 
some neighborhood of the point ( ai , aa , • • • , ap) . By Theorem II, then, the 
functions 2/1 > 2/2 > • • • > yn are linearly dependent in this neighborhood, i. e., 
there exist constants Ci , C2 , • • • , c„ , not all zero, such that the function 

2/ = Ci 2/1 + C2 2/2 + • • ■ + Cnyn 

is zero in the said neighborhood of ( ai , a^, • • • , ap) . But then this func- 
tion vanishes at ( ai , a2 , • • • , a? ) together with its derivatives 2/^^^ , 2/^^^ > ■ • • > 
2/^''^ , and therefore, by hypothesis, vanishes identically. Consequently the 
functions 2/1 > 2/2 , • • • ,yn are linearly dependent in A . 

We now assume the other alternative, viz., that all of the (n — l)-rowed 
determinants of M which contain elements of the first row of M vanish iden- 
tically; in particular, those formed from the first ra — 1 columns of M do. 
We recall that, by hypothesis, if all the I'-rowed determinants formed from 
the first V columns of M and containing elements of the first row of M vanish 
identically, then the determinants formed from the first v columns of ilf', 
and containing elements of the first row and of y — 2 other rows of M , vanish 
identically. We consider then in succession the first n — 1 columns of M , 
the first n — 2 columns, the first w — 3 columns, etc. Suppose that for the 
first n — i columns of M all of the (n — ■t)-rowed determinants which con- 
tain elements of the first row of M vanish identically. We know that this 
happens f or i = and i = 1 . If it happens for i = n — 1 , we shall have 
2/1 = 0, and then yi, y2, • • • , 2/n would be linearly dependent. Suppose then 
that i < n — 1. Then there must exist a particular i,l^i^n — 2, such 
that all the (n — i) -rowed determinants formed from the first n — i columns 
of M and containing elements of the first row of M vanish identically, whereas 

* If for a function of p variables all the partial derivatives o' the first order exist and are 
bounded in the p-dimensional region A , the function is continuous in A with respect to all 
the variables taken together. Cf. Osgood, Lehrbuch der Funktionentheorie, vol. I, second 
edition (1912), page 56. 
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the (n — i — 1 )-rowed determinants formed from the first n — i — 1 columns 
of M and containing elements of the first row of M do not all vanish identically. 
But the identical vanishing of all the (n — i)-rowed determinants of the 
kind considered, formed from the first n — i columns of M , carries with it, 
by hypothesis, the identical vanishing of all the {n — i) -rowed determinants 
formed from the first n — i columns of M' and containing elements of the 
first row and of ft — i — 2 other rows of M . We may therefore apply the 
special case of the present theorem which we have already proved, to deduce 
the linear dependence of 2/1 > 2/2 . • • • > Vnr-i • Consequently, the functions 
2/1)2/2, • • • ,yn are also linearly dependent, and our theorem is proved. 

It is to be noted that in the hypothesis of Theorem IV, a restriction has 
been placed upon the nature of the matrices M and M' which we did not 
find it necessary to impose in the previous general theorems. As may readily 
be seen by an example, it is only in special cases that the identical vanishing 
of all the j'-rowed determinants formed from v columns of the matrix M and 
containing elements of the first row of M implies the identical vanishing of 
all the j'-rowed determinants formed from the corresponding v columns of M' 
and containing elements of the first row and of i' — 2 additional rows of M . 
Probably the most important cases of this type arise when the functions 
2/1 J 2/2 J • ■ • , 2/" ^re solutions of the same completely integrable system of 
partial differential equations. We shall have occasion, in a subsequent sec- 
tion, to apply Theorem IV in just such a case. 

3. Completely integbable systems of homogeneous linear partial 

differential equations 

A system of partial differential equations in one dependent variable is said 
to be completely integrable if there exists one and only one integral of the 
system satisfying arbitrarily chosen initial conditions involving a finite number 
of arbitrary constants and arbitrary functions. The completely integrable 
systems which we shall consider are such that a general solution is expressible 
linearly, with constant coefficients, in terms of a fundamental system of solu- 
tions. Thus, these completely integrable systems, which we shall suppose 
to contain a single dependent variable and any number of independent vari- 
ables, form a natural generalization of the ordinary linear homogeneous differ- 
ential equation. 

So far as we are aware, the only completely integrable systems containing 
partial differential equations of order higher than the first which have hitherto 
been studied have analytic coeflicients and solutions.* From the function- 
theoretic point of view these are of course the most important completely 

* Cf., for example, C. Riquier, Les systemes d'equations aux derivees partielles, Paris, Gauthier- 
Villars, 1910. 
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integrable systems. They have found a place not only in analysis, but in 
geometry as well. The Gauss-Codazzi equations which form the basis for 
the metric differential geometry of surfaces constitute a completely integrable 
system, and in a later section we shall see how important are such systems in 
projective differential geometry. 

The simplicity of the analytic case, if we disregard questions of convergence, 
is apparent. Let us consider, as always, a completely integrable system of 
linear homogeneous partial differential equations containing a single de- 
pendent variable; then by successive differentiation of the given equations, 
all the derivatives of the dependent variable are expressible uniquely as 
linear combinations of a finite number n among them, so that by well-known 
methods it is seen that the general solution contains exactly the same number 
n of arbitrary constants, which enter linearly. 

We shall confine our attention throughout to linear homogeneous com- 
pletely integrable systems in one dependent variable, though most of what 
we shall have to say may readily be extended to the case of several dependent 
variables. We shall not, however, suppose the coefficients in our equations 
to be analytic* 

Let us first construct a completely integrable system which shall have as a 
fundamental system of solutions an assigned set of functions. Suppose that 
the n functions 2/1,2/2, • • • , ?/« of the p independent variables ui, U2, • • • , Up 
possess enough partial derivatives, of any kind or order whatever, to form a 
determinant of the nth order. 



(3) 



W^ 



yi, 


2/2, 


} 


yn 


y'P, 


y'P, 


y 


2/i'' 


yf\ 


yf. 


J 


2/L^' 


2/r^ 


ifr'\ 


> 


y':-'' 



which vanishes nowhere in a region A , and in which, as heretofore, we have 
denoted derivatives of a function y by y^^^ , y^^^ , etc., the indices being no 
indication of the kind or order of the derivative. Suppose also that all of 
the derivatives of the first order exist for all of the elements in the determinant 



* Completely integrable systems in which the coefficients are not restricted to being analytic 
were curiously enough the first to be studied, but only such systems as consist of equations 
of the first order. These systems are really systems of total differential equations, and the 
classic work of Mayer in connection therewith is well known. See his memoir, Ueher un- 
heschranki integrable Systeme von linearen totalen Differentialgleichungen und die simultane 
Integration linearer partieller Differentialgleichungen, Mathematische AnHalen, 
vol. 5 (1872), pp. 448-470. 
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W . Then we may set up as follows a system of partial differential equations, 
of which each of the functions yi,yi, • • • , ?/„ is a solution, and such that any 
solution whatever, say ^ , is of the form 

^ = Ci 2/1 + C2 t/2 + • • • + CnVn, 

where the c's are constants. Let us write down the system of n'p equations 
dy 



dUk 



= dS' *> 2/ + af' *' 2/(" + d!' "'' 2/(2' + • • • + a^^' 2/("-« , 



^ = oo. *) 2/ + o?. *) 2/(1) + 4' "-> 2/(« + • • • + dl-^i y^^-'-> , 



(4) ^y^ 

duk 



= af *> 2/ + di' *> 2/^" + af *' 2/^'^ + • • • + «»-? 2/^""" • 



32/' 



(«-i) 



Swfc 



= «(»-'• *'y + ol"-''*'2/^i) + 4"-'-*) 2/^2) + . . . + a<rri''«2/*""". 



(fc = l,2, •••, p). 

We may so determine the coefficients di' *' of this system that the functions 
2/1,2/2) • • • » 2/» ^ill be solutions of the system. If in any one of the equations 
these functions be substituted successively, n linear equations are obtained 
which may be solved for the coefficients a of the first mentioned equation 
since the determinant W is different from zero. If the coefficients a are deter- 
mined in this way, the system of differential equations will surely have as 
solutions each of the functions yi, y^, • • • , yn- In fact, the system of np 
differential equations may be written in the compact form 



(5) 



2/1. 



2/2, 



2/ni 



2/ 



2/i'\ 


yf, ■■ 


> yn > 


2/(1) 


2/f, 


y'i\ •■ 


if 

> yn > 


2/(2) 


2/^^ 


yf-'\ •■ 


-, yr\ 


y^n-1 


dy? 


dy'^ 


djM' 


dy^'^ 


duk ' 


duk' 


' duk' 


duk 






(J=0 


,1, •• 



= 



n-l; k = 1,2, ••• , p) , 

which makes the statement evident. Moreover, it is also evident that any 
expression of the form 



y = Ciyi + CiPi + ••• + Cnyn 



1916] 



LINEAE DEPENDENCE OF FUNCTIONS 



497 



will also satisfy the system. Lastly, any solution of the system must be of 
this form. For, let y be any solution of the system, and consider the two 
matrices 

M = Mn-i{yi, y-i, •■■ ,yn,y), 
(6) 

M' = Mg{yi,y2, •■• ,yn,y), 

where M' differs from M by the additional rows consisting of all of the partial 
derivatives of the first order of the elements of M which do not already appear 
in M . The derivatives appearing in these additional rows are precisely those 
which appear in the left-hand members of the system of differential equations 
(4), so that each of the additional rows is a linear combination of the n rows 
of the matrix M . Consequently, all of the ( n + 1 )-rowed determinants of 
M' which contain the elements of M vanish identically. But one of the n- 
rowed determinants of M , namely W , vanishes nowhere in ^ , so that by 
Theorem II the n ■\- I functions 2/1,1/2, • • • , yn, y are linearly dependent, 
and in fact 

y = Ci 2/1 + C2 2/2 + • • • + C„ 2/n • 

We may state our result as follows: 

Theorem V. If the system of partial differential equations (4) has n solutions 
yi, y2, • • • , yn for which the determinant (3) vanishes nowhere in A , then any 
solution of the system will be linearly dependent on these n solutions, and any 
function linearly dependent on these n solutions will also be a solution. 

We have therefore obtained a completely integrable system of which the 
given functions 2/1,2/2, • • • , 2/n constitute a fundamental system of solutions. 
The form in which we have chosen to write the system of differential equa- 
tions is not generally the simplest for any given particular case; in fact, in 
many cases most of the equations will be found to be merely identities, and 
others of them obtainable from the rest by mere differentiation. The com- 
pletely integrable system will in such a case reduce to fewer than np equations. 
An example will make the point clearer. Suppose that the four functions 



yk = yk{ui, U2) 



{h = 1,2,3,4) 



are such that the determinant 



W = 
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du^' 
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is zero at no point of the region A . Then the completely integrable system 
of which the four functions yi, y^, yz, yt form a fundamental system of solu- 
tions will be 

dy dy dy _ dy d^ y _ d^ y 

dui dui ' dUi du2 ' dui du^ dui du2 ' 

d^y , i^y , ^y , J ^^y 

aui dui dUi dui dui, 

d"^ y , , , dy , dy „ 9^ w 

5-1 = a' y + b'^ + c' ^ + d' ^-^ , 
out dui du2 aui dUi 

^'y =a"y + h"^+ c"^+ d" ^'^ 
dui du2 dui du2 dui dU2 ' 

JlJL^ = a'" y + h"'^ + c'"^ + d'" -^^^ ■ 
dui dui dui du2 dui du^ ' 

we have of course presupposed the existence of all the derivatives which we 
have written down. The first three of the above equations are identities, and 
as a matter of fact the third of them ought to occur twice. The completely 
integrable system therefore consists essentially of four instead of 2 X 4 = 8 
equations. Moreover, if dd' — 1 is nowhere zero in A , the last two of the 
differential equations may be obtained from the two immediately preceding 
them by differentiation and algebraic processes. The completely integrable 
system of eight equations in this case reduces therefore to the two equations 



d'^y ,K^y_L_ ^y , 1 ^"^y 

-r-2 = ay + 0^ — + ^a — + "5 5" ' 

dui dui dU2 dui 0M2 



VT = a' y + V -r— -\- c' -r— + d' T — 5— . 
dui dui du2 duidu2 

It is evident that in any particular problem this form of the system is most 
useful; however, in the very general questions which we shall treat subse- 
quently, the form (4) in which we first wrote the system seems of great assist- 
ance, in spite of the fact that many of the equations may be superfluous. 

From what we have already seen, the analogy of homogeneous linear com- 
pletely integrable systems to ordinary homogeneous linear differential equa- 
tions in very close. In later sections the similarity will be made to appear 
even more striking; at present we take occasion to make the following remark. 
If in system (4) we keep k fixed, we have a system of n differential equations 
which may be regarded as a system of ordinary differential equations in which 
the independent variable is Uk and the dependent variables are n in number, 
viz., y , y'-^^ , 2/'^' , • • • , j/^""^^ . The completely integrable system (4) may there- 
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fore he regarded in p ways as a system of n ordinary differential equations in n 
dependent variables, containing p — 1 parameters. 

4. The integrability conditions and wronskian of a completely inte- 

grable system 

In the last section, the completely integrable system (4) was set up on the 
supposition that a fundamental system of solutions was known. Naturally, 
if a system of this form, which we shall write 

(7) -f— = ZaV'*'2/''^ U = 0,l,---,n-l;k = l,2,---,p), 

OUk i=0 

is to be completely integrable and have precisely n fundamental solutions, 
the coefHcients cannot be arbitrary. We shall in the present section investi- 
gate certain relations which exist among the coefficients of system (7), under 
the supposition that the system has a fundamental system of n solutions 
2/1 > 2/2 , • • • , 2/„ in the region A . 

Let us call the quantities y , y^^^ , y^''^ , • • • , t/'""^^ which occur in the right- 
hand members of equations (7), the primary derivatives of the function y. 

The determinant PF given by equation (3) is formed of the primary deriva- 
tives of a fundamental system of solutions of the differential equations. We 
shall call this determinant the wronskian of the system of solutions. In 
setting up the completely integrable system (4) we assumed that the wronskian 
of the fundamental system of solutions vanished nowhere in A ; this is equiva- 
lent to the assumption that if for every solution of the system of differential 
equations a relation exists among the primary derivatives, of the form 

(8) aoy + ai y(» + a^ 2/^^' + • • • + a^i ^f"-" = , 

then all of the functions olo, ai, • • • , a„_i must he identically zero. 

We shall from now on assume that in the region A the coefficients aV' *' of 
system (7) are bounded and possess all partial derivatives of the first order. 
Then since by supposition there exists a function y not identically zero which 
satisfies equations (7), we may differentiate these equations and obtain in 
two ways an expression for d^ y''^'>/duk du, {k ^ I) . In fact. 



fc oui i=o \ oui " dui J 



duk du, 

dui dUk i=o V dUk ' dUk J ' 

In the right-hand members of each of these, the derivatives of y^^ may be 
replaced by their expressions in terms of y , y^^^ , • • • , j/^""^^ taken frofn sys- 
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tern (7). In this way two different expressions of the form (8) are obtained 
for a derivative d"^ y^^'^ /dukdui (k 4= I)- But these two expressions are to 
be identically equal, so that by our assumption the coefficients of y^'^ in the two 
expressions must be identically equal. We may therefore state the theorem: 
Theorem VI. A necessary condition that system (7) have a fundamental 
system of n solutions is that the coefficients satisfy identically the relations 



(9) 



dui 



+ E af '^ af '^ = ^ + E aV- " di' "' 
( v,j = 0,1, ■•• ,n - 1; k,l 



1,2, 



P). 



We shall speak of equations (9) as the conditions of complete integr ability, 
or as the integrability conditions, of system (7). 

From the set of integrability conditions (9) let us pick out n equations, 
by keeping k and I fixed and choosing v=j = 0,l,--',n — 1. The equa- 
tions are 

-.^- + i:(^'''a!f' = -^ + Eay-^a^;-*'. (j = 0, 1, ••., n - D- 

OUi i=o OUk i=0 

Adding them, we find without difficulty that all terms drop out except the 
first derivatives, and in fact 



that is, 
(10) 



\i daf'''> ''yidaf''> 
y=o dui j=a dUk 

r\ n-l -) ti-1 

^Eay'*> = /-Eay'^ 

duij^o OUk j=o 



(k,l = l,2, ■■■,p). 



Consequently, the n quantities E"=o «?' *' are the first derivatives of some 
function f{ui,U2, • • • , Up) , and we may write 



(11) 



dUk j=o ' 



{k = l,2, ■■■,v)- 



We shall presently find relations (10) and (11) very useful. 

Let us now form the wronskian (3) of any fundamental system of n solu- 
tions 2/1,2/2, •■• ,ynoi system (7) : 



W = 



2/1. 



,,(1-1) 



2/2, 



W2 , 



y^:' 



y':-'' 



We differentiate W with respect to Uk', if the result be written as the sum of 
the n determinants formed from W by replacing a row theteof by the deriva- 
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tives of that row with respect to ujc , and if use be then made of equations (7) 
to replace any derivative dy'^P/duk by an expression linear in yi , y^P , • • • , 
2/'"~'\ it will be found without difficulty that 

^^W = T.af'-> = -^W {k = l,2,...,p). 

Consequently, 

W = const, e' . 

We may therefore state the theorem, which is a direct generalization of the 
familiar theorem of Abel in the theory of an ordinary homogeneous linear 
differential equation of the nth order: 

Theorem VII. The wronskian W of a fundamental system of n solutions 
of the completely integrable system (7) may be determined by a quadrature from 
the coefficients of the system, and is given by the expression 

(12) W = const, ef, 

where f may be found by a quadrature in virtue of equations (10) or (11). 

As a corollary we have the theorem: 

Theorem VIII. If the wronskian vanish at any point of the region A , it 
vanishes identically in A.* 

We remark in passing that equations (10) and (11) might have been deduced 
directly from a consideration of the wronskian. For, equations (10) may be 
deduced by differentiation of the relations 



(13) ^^WT.af^ (k = l,2,---,p). 

duk j=o 



But the relations (13) are of interest in themselves, because they show that 
the function S^=o a"'*^ is the derivative, with respect to ma;, of a function/, 
ioT k = 1 , 2 , • • • , p , and this without necessitating the assumption of the 
existence of any of the first derivatives of the coefficients di' *' of system (7). 

5. The normal form op the completely integrable system 

A necessary condition that a system of partial differential equations of 
the form (7) have a solution not identically zero is that the coefficients of the 
system satisfy the integrability conditions (9). Whether the condition is 
also sufficient remains to be seen; as a matter of fact, however, an example 
will show that a system of form (7), whose coefficients satisfy the integrability 
conditions, need not have a solution which is not identically zero. For, if 
we select as primary derivatives y and d^ y/du^ , and suppose the independent 
variables to be u and v , we may write down the system of form (7) : 

* We denote by A the region in which all of the solutions 2/1,2/2, • • • , Vn exist. 
Trans. Am. Math. Soc. 33 



502 G. M. 6EEEN: [October 

% , .^'^y dy , , 7,9^ y 

a„3-« y + t> Q^2, Q^2Q^ « 2/ + » a,,2. 

which has no solution different from zero even if the integrabiHty conditions 
of form (9) are satisfied for the system, unless of course other relations exist 
among the coefficients. In fact, differentiation of the first equation of the 
system twice with respect to u will lead to two new equations which are in 
general incompatible with the first and third equations of the system. In 
our example this results from the peculiar choice of primary derivatives, a 
gap occurring between y and 3^ y/du^ . 

We shall presently show that it is always possible, for a system of form (7) 
which is known to have n fundamental solutions, to makie a particular choice 
of primary derivatives which removes the above difficulty, although at a 
slight sacrifice in generality. Let us first define what we shall call a chain of 
derivatives and a normal set of derivatives. A chain of derivatives of a func- 
tion 2/ is a set of derivatives y , y^^'> , 2/'^' , • • • , 2/'''' , such that y'-'^ is of the ith 
order and may be obtained by a differentiation from y^^"^^ . Thus, denoting 
partial differentiation by subscripts, we see that the quantities y, yu, yuv, 
yuuv form a chain, while y ,yu, yuv, 2/««« do not because yuuu cannot be formed 
from yuv by differentiation. 

A normal set of derivatives is a set such that each element of the set belongs 
to at least one chain which can be formed from elements of the set. The 
normal set need not consist of a single chain, and a given member of the set 
may belong to more than one chain. 

We shall now show that if system (7) has exactly n fundamental solutions 
2/1)2/2, • • • , yn, the primary derivatives may under certain conditions be 
so chosen as to form a normal set. However, in order not to interrupt our 
reasoning later, we shall first prove a lemma which is to be used in the course 
of the discussion. 

Lemma. Suppose that system (7) has n solutions 2/1,2/2, • • • , yn whose wron- 
skian W vanishes nowhere in the region A . Then 2/1,2/2) ' ' ' > yn cannot all be 
solutions of a system R of form (7) having less than n primary derivatives, these 
primary derivatives being selected from among those of system (7). 

Let y , 2/'^' , y^^^ , • • • , 2/'*^' (r<7i — l)be the primary derivatives of the 
system R . Then the lemma follows at once if in the matrix 

M = Mr{yi,y2, ■ ■ ■ , Vn) 

at least one of the (r + 1) -rowed determinants vanishes at no point of A, 
for then this determinant plays the role of a wronskian for system R . Sup- 
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pose, however, that each of the (r + l)-rowed determinants vanishes some- 
where in A. They can not all be identically zero in A, because then the 
wronskian of 2/1 , 2/2 , • • • , Vn would be zero also, contrary to hypothesis. 
Consequently one of these determinants, say that formed from the first r + 1 
columns of M , is different from zero at some point of A , viz., ( ai , • • • , a^,) , 
and therefore, since it is continuous, it differs from zero in a p-dimensional 
region A' containing this point. Then in this region A' , every solution of 
system R is linearly dependent on yi, 1/2, • • • , t/r+i , so that in A' every solu- 
tion of (7), being also a solution of R, will be hnearly dependent on yi, 2/2, 
• • • > 2/r+i • Then in particular we shall have in A' 

yr+2 = Ci 2/1 + C2 2/2 + • • • + Cr+1 yr+1 , 

the c's being constants not all zero. Therefore the function 

y = yr+2 - Ci 2/1 — C2 2/2 — • • • - Cr+1 yr+1 

will vanish identically in A' , and so will its derivatives y'-^^ , y'-^'^ , • •• , y'-"^ • 
But it is easily seen that if system (7) has n solutions 2/1 > 2/2 , • • • , ?/« for which 
the wronskian vanishes nowhere in A , then any solution all of whose primary 
derivatives vanish at some point of A will vanish identically in A . Conse- 
quently the function y considered above vanishes identically not only in A' , 
but throughout A . Therefore we conclude that system (7) would have only 
r + 1 linearly independent solutions, which contradicts the hypothesis. 

We now proceed to show how system (7) may be reduced to a similar 
system for which the primary derivatives form a normal set, under the hypo- 
thesis that system (7) has exactly n fundamental solutions yi, yi, ■ • • , yn 
whose wronskian W does not vanish at any point of a region Aoip dimensions. 
We assume, then, that the n primary derivatives of system (7) do not form a 
normal set. Then we may select from among these primary derivatives a 
normal set which consists of at least one element, namely y; let this normal 
set be 2/ , y^^^ , 2/^^^ , • • • . 2/'''' , {r < n — I) . Suppose, furthermore, that it is 
the largest normal set which may be formed from among the primary deriva- 
tives. Then of the first derivatives of y , y'-^^ , y^^^ , • • • , y^^^ none is a primary 
derivative of system (7), for otherwise our normal set could be enlarged.* 
Now, all of the first derivatives just mentioned are expressible linearly in 
terms of the primary derivatives of system (7), the coefficients in the expres- 
sions being continuous in A . This gives us p ( r -)- 1 ) equations, in the right- 
hand members of which appear some or all of the quantities y , 2/*^^ , 2/*^' . • • • j 

* Lest a misunderstanding occur here, we point out that if j/"* and j/'"* are primary deriva- 
tives, and if dy^''>/duj = j/'"', then this equation is one of the equations of system (7); also, 
if 2/*" be one of the normal set y, 2/<w , • • • , j/''' under consideration, then y^'"'> belongs to the 
said normal set also. 
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y^^^ . But in the right-hand member of at least one of these equations must 
appear some one of the primary derivatives, say 2/''^ , which is not one of 
our normal set, for otherwise the p{r + 1) equations would constitute a 
system of form (7) whose primary derivatives are y , y^'^^ , y^^^ , • • • , y^^^ , i. e., 
the said p{r + 1) equations would constitute a system R such as is con- 
sidered in the lemma proved above, and the original system (7) could then 
have no more than r + 1 linearly independent solutions, contrary to hypothesis. 
Consequently, of the first derivatives of y , y^^^ , y^^^ , • • • , y^^^ , at least one 
has a linear expression, in terms of the primary derivatives, which contains a 
primary derivative y^^^ which is not a member of the normal set y , z/'^' , 2/'^' , 
• • • , 2/^''' • In other words, the coeflBcient of y'-"^ in this expression is not 
identically zero; since this coefficient is continuous, there must exist a p- 
dimensional sub-region of A , say Ai, at no point of which it vanishes. Con- 
sequently, in the region Ai we may express the primary derivative t/^'^ linearly 
in terms of all the other primary derivatives, and the additional quantity 
which occurs in the left-hand member of the equation from which we have 
just found the expression for 2/''' • This additional quantity, we recall, is 
one of the first derivatives of one of the elements of the normal set y , y'-^^ , 
y(.i) , . • . , y^^> , and is not itself a primary derivative. Let us denote it by y^'^ . 
Then throughout (7) we may replace y'-'^ everywhere by its equivalent ex- 
pression, and thus obtain a new system of differential equations, of form (7), 
in which the primary derivatives are the same except that 2/^'^ has taken the 
place of 2/''^ • This new set therefore contains in it a normal set composed of 
the previous normal set y , y^^^ , y^^^ , • • • , y'-''^ and the additional derivative 
y'-'^ . The new wronskian of the system of differential equations vanishes 
nowhere in the p-dimensional region Ai. But although the coefficients in 
system (7) are supposed to have all of their first derivatives, it is logically 
possible that all of the first derivatives do not exist for the coefficients of the 
new system. We shall suppose that these derivatives do exist.* Then on 
the new system of differential equations we may proceed as before, enlarging 
step by step the number of primary derivatives which form a normal set until 
finally a system is obtained for which all of the n primary derivatives form a 
normal set. In this reduction we have supposed of course that certain deriva- 
tives exist for the coefficients in system (7). In any event, it is clear that 
such derivatives need not exist of order higher than n — r — 1 , where r is 
the number which appears in the above discussion. 

We shall henceforth suppose that the primary derivatives in system (7) 



* The writer does not know whether they exist as a consequence of the previous hypotheses 
on system (7) or not. Since the form to which we shall ultimately reduce our system seems 
essential for the proof of the existence theorem, it would be of interest to have this point 
settled. 
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form a normal set, and shall then speak of system (7) as being in the normal 
form. The chief importance of the normal form lies in the fact that system 
(7) in that form may be regarded as a system of np partial differential equa- 
tions of the first order in the p independent variables ui, u%, ■ • • , Up, and the 
n dependent variables y , y^^^ , z/'^' , • • • , y'-"~^'> , since among equations (7) 
occur a certain number which define each primary derivative as the first 
derivative of another primary derivative belonging to the same chain.* 

6. The canonical form of a completely integrable system 

It is well known that by a transformation of the dependent variable of the 
form y = 'K{x)y any ordinary homogeneous linear differential equation of 
the nth order, 

d" y d"~^ y 

whose coefficients are continuous functions of x , may be transformed into one 
of the same form in which the coefficient pi is zero. We shall now show that 
a similar thing may be done for certain completely integrable systems. 

We introduce a new dependent variable, y , connected with the old by the 
equation 

y = >^y> 

where X is a function oi ui,U2, • • • ,Up. Any derivative of y , say y'-^^ , is con- 
nected with derivatives of the new dependent variable, ^ , by an equation of 
the form 

where the part omitted in the right-hand member is linear in derivatives of y , 
and contains all those derivatives of y from which y'-^^ may be obtained by 
differentiation. 

Let us now carry out this transformation on a particular kind of completely 
integrable system, namely, one in which the primary derivatives not only 
form a normal set, but have the additional property that, if a derivative y'-'"' 
belong to the set, then every derivative of lower order from which y'-^'> is ob- 
tainable by successive differentiation is likewise a primary derivative. It is 
not difficult to see that in this case, if y be replaced everywhere by X^ , the 
system of differential equations is changed into a new system, with the de- 
pendent variable y , which new system may be thrown into exactly the same 
form as the original system. Let us write the new system 

(14) f^=Zai^'*'^« 

(i = 0, 1, ■■■ ,n -I; k = l,2y---,p), 



* See footnote on page 503. 
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where the coefBcients ai-''*^ are functions of the old coefficients df''^ (of 
system (7)) and also of X and its successive derivatives. Of course the actual 
expressions depend upon the form of each individual primary derivative, 
which form we have not specified. 

We may, however, obtain explicit results with regard to the transformation 
of the wronskian and its first derivatives. Let us denote by W the wronskian 
for the transformed system (14). We recall that the set of primary derivatives 
has the property that if y'^^^ belong to the set all the derivatives of y from 
which 2/^^^ may be obtained by differentiation also belong to the set. Conse- 
quently the wronskian W for the new system is easily seen to be related to the 
wronskian for the old system by the equation 

(15) W = \^W . 

Differentiating this equation with respect to Uh , we have 

dw ^ dw , , — ax 

auk duk auk 

or, using equations (13) and (11), 

aUk oUk aUk 

where / denotes the function which takes the place of / for the new system of 
differential equations (14). Substituting for W its value as obtained from 
(15), and recalling that W is nowhere zero in the region A , we find that 

\f = \^-n^ (fc = l,2,...,p). 

dUk OUk OUk 

This result enables us to state the theorem : 

Theoeem IX. Suppose a system of partial differential equations of the 
form (7) has the following properties : 

1°. It is completely integrable, with n fundamental solutions whose wronskian 
(3) vanishes at no point of A . 

2°. The set of primary derivatives is such that, if 7/^^' be any one of the set, 
then all the derivatives of lower order from which y'-^^ may be obtained by differ- 
entiation also belong to the set. 

3°. All of the primary derivatives exist for each of the np coefficients af'^^ 
(i = 0,l, ■■■,n-\; k = 1,2, ■■■,p). 

Then the system of differential equations (7) may be transformed into another 
uniquely determined completely integrable system (14) /or which all of the quantities 

~- =llaf"'' (fc=l,2, •••,p), 

OUk >=o ■' 
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are zero, by the transformation of the dependent variable y = \y , where 

(16) X^ - nx~ =0, X = const, e^j'", 

aUk aUk 

The system (14) thus obtained will be unique, and its coefficients a^/'*' 
will be explicit expressions in the coefficients of system (7) and their deriva- 
tives. Let us call this uniquely determined system the canonical form of 
system (7). 

It should be noted that part 3° of the hypothesis is necessary in order that 
X^ may be substituted for y throughout equation (7). 

Let us return to the general transformation y = \y , which converts system 
(7) into system (14); we now suppose that system (14) is not in the canonical 
form. The coefficients of system (14) are expressible explicitly in terms of 
the coefficients of (7) and their derivatives, and X and its derivatives. Con- 
sider any function F of the coefficients of (14) and their derivatives. Let F 
be the same function of the corresponding coefficients of (7) and their deriva- 
tives. The function F is equal identically to a function of X and the coefficients 
of (7), and the derivatives of these quantities. If in particular this latter 
function is related to F by the identity 

F = VF, 

where r may or may not be zero, we call the function F a seminvariant of 
system (7). 

The following theorem is a consequence of the fact that system (7) and the 
transformed system (14) have exactly the same canonical form : 

Theorem X. The coefficients of the canonical form are seminvariants of 
system (7). All other seminvariants are functions of these and of their derivatives, 
if such derivatives exist. 

The reduction of a completely integrable system to the canonical form is 
seen by the above discussion to be a purely mechanical process, performable 
in any case coming under the hypotheses of Theorem IX. This fact is of 
practical importance in projective differential geometry, if the general method 
of Wilczynski be followed. We shall return to this question in the next sec- 
tion, applying the method to a particular geometric problem. 

In § 8, we shall give an existence theorem for systems of form (7) in which 
the primary derivatives constitute a normal set. It should be noted that 
in the hypothesis 2° of Theorem IX, however, the set of primary derivatives 
is still further restricted; as a matter of fact, we may show by an example 
that this restriction is in general necessary, so that not every completely 
integrable system may be reduced to a unique canonical form. Denoting 
by subscripts partial differentiation with respect to the two independent 
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variables u and v , let the primary derivatives he y , yu, ym ■ These form a 
normal set. If yi, y^, ys be functions oi u, v for which in the region A the 
wronskian 

W = \y, yu,yuv\ 

is nowhere zero, we may set up the completely integrable system having 
yi> 1/2, 2/3 as a fundamental set of solutions: 

y„ = ay + byu + cyuv , 

yuu = a' y + b' y„ + c' yuv , 

yuuv = a" y + h" y^ + c" yuv , 

ymv = a'" y + V" yu + c'" yuv , 

where in general c 4= . Making the transformation y = 'Ky , we find that 
W is transformed into 

W = |X^, \yu + \uy, X^„» 4- X» yu + Kyv + X^ y\ = (X^ + cX^X„) W , 

in which use has been made of the first of the system of differential equations. 
Now, if c 4= — a condition, which may be realized at pleasure — this relation 
between W and W is not of the form (15), so that the differential equation (16) 
which determines X is replaced by one of the second order. Therefore, 
although the system may be reduced to the canonical form, this form will not 
be unique. 

Of course, instead of making the said hypothesis 2° of Theorem IX, we 
might demand merely that the wronskian and its transform be related by 
equation (15). Our example would then be included in this case, if c were 
identically zero in A . 

7. Geometbic application; curvilineae coordinates* 

Completely integrable systems of the kind we have been considering have 
attained increased importance of late owing to their systematic use in certain 
fields of projective differential geometry. Since any solution of the system is 
linearly dependent on any fundamental set of solutions, it follows that a geo- 
metric configuration defined by a fundamental set of solutions is a projective 
transformation of the configuration defined by any other fundamental set of 
solutions. Consequently, any geometric property expressed in terms of the 
completely integrable system, being common to all the projective trans- 
formations of the configuration defined by a fundamental set of solutions of 
the system, is a projective property. Wilczynski has developed a general 
method for dealing with certain questions from this point of view, and it 

* This section is independent of § 8. 
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is in connection with this method that the results of § 6 find an important 
application. 

Let us consider in particular the general subject of curvilinear coordinates 
in space of n dimensions.* If j/i , 2/2 , • • • , y.n+i be interpreted as homogeneous 
coordinates of a point in this space, the equations 



(17) 



Vi = ViiUl, U2, ■■■ , Un) 



(1 = 1,2, ...,n + l) 



will define a system of curvilinear coordinates, provided the ratios of the r/'s 
do not all reduce to functions of fewer than n of the independent variables 
Ui, Ui, • • • , Un. We shall suppose that all derivatives of any order exist for 
the functions yt . Since at least one of these functions, say yn+i , is different 
from zero at some point, it will be different from zero in an n-dimensional 
region about this point. We wish, then, the analytic condition that the n 
ratios rji = yi/y„+i, {i = 1 , 2 , • • • , n) actually define 00" points. This 
condition is the non-vanishing of the jacobian of the functions rji, and this 
may be reduced to the non-vanishing of the determinant 



(18) 



W=' 



yu 



2/2, 



2/»+i 



y'P, 


yV; ■■ 


•, 2/i'i. 


y?\ 


2/'/', •• 


) 2/«+i 


y'f, 


2/?', •• 


•, y':lx 



where y^^ , 'ifi' , etc., denote dyi/duj, dyi/duj, etc. We may therefore set 
up a completely integrable system of form (7) of which the n + 1 functions 2/i 
constitute a fundamental set of solutions, and in which the primary deriva- 

(n) 



tives are y , 2/'^' , 2/^^* 



(19) 



dy 

duk 



yv 



= /"': 



This system is 

^ = Eal^-'*>2/^^> 



U,k = l,2, 



,n). 



Let us differentiate W with respect to Uk', then on using equations (19) we 
find without difficulty that 



dW 
dUk 

We may therefore put 



= (a\''*' + dr"' + •■■ +<"■*') ^ 






It 
dUk 



{k =1,2, ••-,«). 



* The case n = 2 , i. e., nets of plane curves, has been studied by Wilczynski, One-parameter 
families and nets of plane curves, these Transactions, vol. 12 (1911), pp. 473-510. 
The present writer considered the case n = 3 in his Columbia dissertation. Projective differ- 
ential geometry of triple systems of surfaces, Lancaster, Pa., Press of The New Era Printing 
Company, 1913. 
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Now, there is an arbitrary feature introduced in equations (17) because of 
the homogeneity of the coordinates. If X be any function oi ui, u^, • • ■ , u„ , 
then the points defined by the coordinates yt (i = 1,2, • • ■ , n + 1) will be 
the same as those defined by equations (17), if yt = Xj^i . Consequently, if a 
projective property of the configuration is to be expressed in terms of the 
coefficients of the completely integrable system (19), it must be in terms of 
what we called in § 6 the seminvariants of the system, viz., those combinations 
of the coefficients of equations (19) and their derivatives which remain essen- 
tially unchanged when the system is subjected to an arbitrary transformation 
of the form 

y = ^y- 

Let us carry out this transformation of the dependent variable, thus con- 
verting system (19) into one of the same form, 

(20) ^^f'K ^-=Za?'*^^W ij,k = l,2,...,n). 

The coefficients of this system are functions of the coefficients of (19) and 
their derivatives, and also of X and its derivatives. In particular, 

af "> = a'/' "^ - X<*> A (i + ^ ) > af '^ = df' ^^ - 2X(*>/X 

U, fc = 1,2, •••, n), 

so that, on denoting by df/dut the quantity S"=i 5^'*^ we find directly that 

auk ouk aUk 

This also follows from the general results of § 6. Consequently, the quantities 
J(i) = df/duk will all be zero if X be chosen as a solution of 

Xf*=) - (n + l)X(*^> =0 (fc = 1,2, ••.,«), 

which gives 

X = const. e-^/(«+i) . 

With this choice of X the coefficients of the new system (20) have the values 






+ l\t[' ^ duj n + 1^ ^ )' 



n + 



ar 



1 / ^ <9f f J'> 1 \ 

„u, « 1 — ± — ( V gfj, K) ni) _ "fl — nj) f (*) 1 



50-,*) = ««.*)_/(.■) (y =^fc; ,-,y,fc = 1,2, ...,n), 
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where throughout single superscripts denote differentiation. The quantities 

(21) are the fundamental seminvariants of system (19), and any semin variant 
of the system is a function of these and their derivatives. 

There still remains, however, an arbitrary feature in equations (17). This 
is in the choice of parameters. Since the system of curvilinear coordinates 
remains unchanged under the transformation 

(22) Uk=Uk{uk), (k = 1,2, ■■■,n), 

it follows that projective properties of the configuration are expressible in 
terms of those seminvariants which remain essentially invariant when system 
(19) is subjected to the transformation (22). It is possible to calculate 
explicitly, and without any difficulty, these invariants, but our object here is 
not to give a systematic theory of curvilinear coordinates. 

The fact that even in this very general problem the invariants may be 
calculated explicitly, makes it highly probable that the analytic methods of 
Wilczynski are practicable in handling problems in the projective differential 
geometry, not only of configurations in spaces of a few dimensions, but also 
in space of n dimensions. The only methods which have hitherto been used 
systematically in studying configurations in hyperspace are those of Segre.* 
The ideal method seems to the present writer to be one in which all the pro- 
jective differential properties are expressible by means of a unique analytic 
apparatus, and this is to be found in the completely integrable system of 
which a fundamental system of solutions defines the particular geometric 
configuration under discussion. 

8. The existence theorem fob completely integrable systems. 
The linear dependence op solutions 

The theorems of previous sections concerning the system of partial differ- 
ential equations (7) presuppose the existence of a fundamental system of n 
solutions. The question of this existence presents no essentially new diffi- 
culty, for the remark made at the end of § 5 renders immediately applicable 
to system (7), which is supposed to be in the normal form, a well-known exist- 
ence theorem concerning systems of total differential equations. f 

* For references to some of this work, see E. Bompiani, Recenti progressi nella geometria 
proiettiva differenziale degli iperspazi, Proceedings of the fifth international congress of mathe- 
maticians, Cambridge, 1912, vol. II, pp. 22-27. 

In this connection should be mentioned also Prof. Wilczynski's Cleveland address, Some 
general aspects of modern geometry, Bulletin of the American Mathemati- 
cal Society, vol. 19 (1913-14), pp. 331-342. 

t Cf., for example, Ch.-J. de la Vallee Poussin, Cours d' Analyse Infinilesimale, 2me ed., 
vol. 2, chap. VII, § 6. The theorem is essentially that of Mayer, given in the memoir cited 
in §3. 
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A simultaneous system of n total differential equations in n dependent 
variables x^ {h = 1 , 2 , • • • , n) and m independent variables yt {i = 1 , 2, 
• • • , m ) is of the form 

(23) dxH^ildpdyi 

(A = 1, 2, •••, w; i = 1,2, ■■■ ,m.), 

the a's being continuous functions of the ar's and y's in a region A . Let us 
write 

^Vk dyk A=i " 5a;s;' 

Then we shall say that system (23) is completely integrable if its coefBcients 
satisfy the identities 

(24) f.f 

5yk ^Vi 

( i , fc = 1 , 2 , • • • , to; A = 1 , 2 , • • • , TO ) . 

We shall now state the existence theorem for this system. 

Theorem. Suppose that the identities (24) subsist for system (23) throughout 
the region A , and let yf-, • • • , y'-^; ^T > ' " " » ^n' ^^ (^ *^* of values of the inde- 
pendent and dependent variables in the Interior of this region. Then there exists 
one and only one set of n functions Xi, • • • , Xn of the m independent variables 
2/1 > ' • ' y Vm, which satisfy the system of differential equations and which respec- 
tively take on the values xf^ , • • • , a;^°' when yi = i/i"' , • • • , ym = 2/m' • 

Consider now system (7), viz., 

duk h'^' y 

(i = 0, 1, •••,«-!; fc = 1,2, ■•-,?), 

in which we suppose that there are exactly n primary derivatives, and that 
they form a normal set. Let us regard y , y^^^ , • • • , y'-"~^'> as dependent 
variables; then the system becomes a very particular case of the system f23) 
which we have just been considering. The integrability conditions (24) take 
precisely the form of equations (9). Moreover, from the remark made at 
the end of § 5, system (7) when regarded as a system of partial differential 
equations of the first order with the dependent variables y , y^^^ , • • • , y'-"~^^ is 
exactly equivalent to system (7) with its original significance, viz., that in 
which y , y^^^ , • • • , j/^""^^ are primary derivatives of the function y , each 
superscript denoting a certain particular, though unspecified, partial deriva- 
tive of y . The exact equivalence of the two points of view follows from the 
fact that the primary derivatives form a normal set. In other words, if 
system (7) be regarded as a system of partial differential equations of the 
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first order in n dependent variables, then a set of solutions y , y^^^ , • • • , y^"-~^^ 
will have the property that y^^^ , y'-^^ , • • • , y^""^' are exactly those derivatives 
of y which the original superscript notation signified; and all this without 
requiring the addition of new equations to system (7). 

The force of these remarks, especially the last, may be seen by means of 
the completely integrable system in two independent variables u and v , with 
primary derivatives y , yu, yv' 

yuu = ayu + hy^ + cy , 

(a) 2/„^ = a' yu + h' y^ + c' y, 

yvv = a" y„ + h" 2/„ + c" y , 

where the subscripts denote differentiation. This system may be replaced 
by a system with three dependent variables by means of the substitution 



ih) 


Vu 


= 2/(1), 


Vv = y^^^ , 


and will then have the form 








2/L'' = 


a2/W +%e^ +C2/, 


(a') 


^f = 


a' 2/« + V 
a' ^») + V 


y^'^ +c'y, 
y^'^ + c'y, 



y':-> = a" y^» + b" y^^^ + c" y . 

This is not in the form (7), however, because all of the first derivatives of 
all of the primary derivatives do not appear in the left-hand members of 
system (a')- But equations (a') and (6) together constitute a system of 
form (7). System (a') alone is not equivalent to system (a), but systems (a') 
and (6') together are equivalent to system (a). 

A better illustration is afforded by the system of form (7) introduced at 
the beginning of § 5, viz., one in which the primary derivatives are y and yuu , 
and which we there saw had no solution different from zero. If we put 

yuv, ^^ y ! 

then the system in form (7) is 

yu = ay + by^^^ , y^ = a' y + b' y^^'> , 

yj) = a"y + b" 2/(1) , 2/(1) = a'" y + V" !/« . 

These four equations, regarded as equations of the first order in the dependent 
variables y , y^^^ , will have sets of solutions y , y^^' different from zero if the 
integrability conditions (9) are satisfied for the system. But the function t/^) 
will not in general be the proper derivative of the function y , viz., 2/m« . This 
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can not happen in the case of the other example, because equations (b), which 
form part of the system in the form (7), make it certain that a solution y , y'-^^ , 
y^^y of the system consisting of equations (a') and (b) will have the property 
that 2/*^' and y^^^ will be the proper derivatives, 2/u and y^ respectively, of y. 
Let us, then, regard system (7) as a system of partial differential equations 



of the first order in the n dependent variables y , y' 



(1) 



y 



(n-l) 



Applying 



the existence theorem for system (23), of which (7) is a particular case, we 
may state the following existence theorem : 

Theoeem XI. Suppose that in the system of partial differential equations 



(7) 



dy 



U) 



dUk 



= Za^'*^2/« 



U = 0,1, ■■■,n-l; k = 1,2, 



P) 



the primary derivatives y , y'-^^ , • • • , y'-"~^^ form a normal set. Suppose further 
that in the region A the coefficients a'-/'*', which are functions of the p inde- 
pendent variables Ui, U2, • • • , Up, satisfy identically the integrability conditions 



(9) 



da\ 



XJ, *) 



dui 



+ Zay'*'a<J''^ = 



dd 



(J, « 



+ ZaF'"at*'*' 



duk 

{p,j = 0,1, ■■■,n -1; k,l = 1,2, ■■■,p). 

, u^p'') be any point of A, and ya, i/o\ • • ■ , 2/0""" be any 
set of n constants. Then there exists one and only one function y of the variables 
ui, U2, • • • , Up which satisfies the system of differential equations, and whose 
primary derivatives y , y^^^ , • • • , j/*""^^ take on respectively the preassigned con- 



Let {uf\ uf>, 



stant values yo , y\ 



,(1) 



, yf ^' at the point ( u'l 



m 



uT, 



<>). 



As in the familiar case of ordinary differential equations, this theorem 
allows us to infer the existence of n linearly independent solutions whose 
wronskian vanishes nowhere in ^ . In fact, we may construct a particular 
set of solutions having this property, by choosing solutions Yi, Y2, • • • , F„ 
such that the primary derivatives of the solution Yi take on the values 

YT = o {j=¥i-i), yv^ = i(i = i-i) 



at the point {u'-i\ 
• • • , F„ is then 



uT, 



u: 



(0) 



) . At this point the wronskian of Fi , F2 . 



1, 

0, 



0, 
1, 



0, 0, 



= 1, 



so that the wronskian vanishes nowhere in /I , by Theorem VIII. Moreover, 
that solution of the system for which the primary derivatives take on the 



values yo , 2/' 



,(1) 



, 2/Sr~" at the point {ui, U2, •••, Up) is given by 
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2/ = 2/0 Ti + 2/0'' Y, + y'PY,+ ■■• + yt~'' F„ . 

If a solution has all of its primary derivatives zero at any point of ^, it 
therefore vanishes identically in A . Let 2/1,2/2, • • • , Vn he any n solutions 
of the completely integrable system; then no function, other than zero, of 
the form 

2/ = Ci 2/1 + C2 2/2 + • • • + CnJ/n, 

can vanish together with its derivatives y'-^^ , y^^^ , • • • , y'-'^~^^ at any point 
of A , unless the constants c, are all zero. This set of functions 2/1,2/2, • • • , yn 
has therefore the properties of the n functions of Theorem III. We may 
therefore apply this theorem to prove the following one, which is also a conse- 
quence of the existence of n solutions whose wronskian vanishes nowhere 
in A: 

Theorem XII. // 2/1 , 2/2 , • • • , yn be any n solutions of the completely 
integrable system of Theorem XI, then the vanishing, anywhere in A, of the 
wronskian of these n functions is a sufficient condition for their linear dependence. 

This theorem is included in the following more general one. They are 
both generalizations of theorems given by Bocher for ordinary differential 
equations.* 

Theorem XIII. /f 2/1 , 2/2 , • ■ ■ , yk {Jc ^ n) be solutions of the completely 
integrable system of Theorem XI, a sufficient condition that they be linearly 
dependent in A is the identical vanishing in A of all those k-rowed determinants 
of the matrix 

M = lf„_i(2/i, 2/2, ••• , 2/t) 

which contain elements of the first row of the matrix. 

The proof of this theorem follows immediately from Theorem IV. The 
matrix M of Theorem IV is in the present case the matrix formed from the 
primary derivatives of the functions 2/1,2/2, ' ' ' , 2/* , while the augmented 
matrix M' has additional rows consisting of the first derivatives of the quanti- 
ties appearing in M . But these first derivatives are linear combinations of 
the primary derivatives, in virtue of the completely integrable system, so 
that if any v columns {v ~ k) he selected from M , the vanishing of all the v- 
rowed determinants formed from these columns and containing elements 
of the first row of M carries with it the vanishing of all the j'-rowed deter- 
minants formed from the same v columns of M' and containing elements of 
the first row of M . Moreover, no solution of the completely integrable sys- 
tem, other than zero, can have all of its primary derivatives vanish simul- 
taneously at any point of A . Consequently all of the hypotheses of Theorem 
IV are satisfied by our functions 2/1,2/2, ■ ■ • , 2/fc , ^^^ the present theorem 
follows at once. 



' M. Bocher, loc. eit., these Transactions, vol. 2 (1901), pp. 139-149, theorem VII. 
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Enough has been given to show the marked analogy between ordinary 
homogeneous linear differential equations of the nth order and the completely 
integrable systems which have come under discussion in this paper. This 
analogy indicates the direction in which facts concerning ordinary differential 
equations may be extended to systems of partial differential equations. It is 
obvious how an adjoint system may be defined; the question naturally presents 
itself as to whether boundary value problems similar to those for ordinary 
differential equations may arise in connection with completely integrable 
systems. Questions concerning analytic functions of several variables suggest 
themselves also. In particular, can the solutions of a given completely 
integrable system be characterized by their properties as analytic functions, 
that is, does there exist a solution of the Riemann problem for such a system ? 
Unfortunately, the theory of functions of several variables has not yet reached 
a degree of development comparable with that of functions of one variable, 
so that it is doubtful whether the generalization is feasible at present. How- 
ever, all of these subjects certainly deserve investigation. 

Harvard University, 
November, 1915 



